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Abstract 

In this paper we evaluate the rcnormalization constants and anomalous dimensions for 
the squark wave function and mass within supersymmetric QCD. These results complement 
the ones obtained in Ref. [Tj and thus provide further confirmation on the applicability of 
dimensional reduction to supersymmetric QCD at three-loop order. The three-loop anoma- 
lous dimension constitute important input to precision predictions of the supersymmetric 
mass spectrum as obtained from the evolution from the GUT to the TeV energy scale. 



1 Introduction 



Supersymmetry (SUSY) (for a review see, e.g., Ref. [2]) has a number of appealing properties 
which classifies it as a promising extension of the Standard Model (SM). Among them are the 
possibility of gauge coupling unification, a dark matter candidate, and a solution to the hierarchy 
problem. 

Although there is yet no clear evidence for the realization of SUSY in nature it is mandatory to be 
prepared both on the experimental and theoretical side. Currently there are several experimental 
groups who eagerly look for signatures of supersymmetry in the data provided by the CERN 
Large Hadron Collider (LHC). As far as theory is concerned it is on the one hand important to 
provide precise predictions for production cross sections involving SUSY particles. On the other 
hand there are a number of quantities which require higher order loop corrections. A prominent 
example is the prediction of the lightest Higgs boson mass which recently became available to 
three loops [3l [H [5] resulting in an uncertainty which can nevertheless be of the order of about 
1 GeV [5]- Another example where higher order corrections within a super symmetric theory are 
very welcome are the renormalization group functions. They are crucial for the running from 
low to high energy scales and constitute an important input for the spectrum generators (see, 
e.g., Refs. [3 [8]) which predict the SUSY spectrum on the basis of only a few assumptions at 
energies of about 10 16 GeV. 

The canonical choice for the regularization used for higher order loop calculations is dimensional 
regularization (DREG). However, it is known since about 30 years that DREG breaks SUSY. As 
a way out dimensional reduction (DRED) has been formulated j9j fT0| [TT] which takes over most 
of the convenient features from DREG and is thus a viable alternative for practical multi-loop 
calculations. It is worth mentioning that DRED is equivalent to DREG for non-SUSY theories 
as has been shown in Refs. [121 [T3l PHI [T5l [TBI [T71 118j . Furthermore it has been demonstrated 
in a number of papers pU EUj EJJ Efl [JJ [23] that DRED is consistent with SUSY QCD at the 
three-loop level. In this paper we provide as new ingredients a further contribution by computing 
three-loop renormalization constants for the mass and mixing angle of squarks in the minimal 
subtraction scheme, which in the context of DRED is called DR. 

The renormalization constants and the corresponding anomalous dimensions up to two-loop 
order has been computed in Ref. |24l [25l [26l [271 [5] . Three- loop corrections have been consid- 
ered in Refs. [191 E01 EU [22], [28] using relations between the beta functions of the gauge and 
Yukawa couplings and the anomalous dimensions of the symmetry breaking parameters that 
can be established in a softly broken supersymmetric theory [291 [30] I31j . In Ref. [JJ the wave 
function renormalization constants of quarks, squarks, gluons, gluinos, ghosts and e scalars and 
the renormalization constants for the quark and gluino mass were calculated to three-loop order 
in the framework of SUSY QCD. In Ref. [JJ also the /3 function for the strong coupling constant 
has been derived from all possible three-point functions. The fact that in each case the same 
expression has been obtained provides a check on the consistency of DRED with gauge invari- 
ance and supersymmetry. In this paper the squark renormalization constants are computed to 
three loops using the component field approach. The main difficulty of this calculation in con- 
trast to the renormalization constants for the gluino and quark masses is that the squark mass 
renormalization constant depends on the masses of the occurring particles in the loops although 
a renormalization scheme based on minimal subtraction is adopted. Furthermore, there is an 
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interplay of the renormalization of the e scalar and the squark mass which will also be discussed 
in this paper. 

The remainder of the paper is organized as follows: In the next Section we derive formulae for 
the squark renormalization constants and briefly outline the procedure used for the construction 
of the exact mass dependence. Furthermore, the renormalization of the e scalars is discussed in 
detail. Our results are presented in Section and Section 2] contains the conclusions. 

2 Formalism 

The calculations in this paper are performed in the framework of SUSY QCD with n q = 5 
massless quarks and a massive top quark (rrit). The scalar super partners of the latter has two 
mass eigenstates (m^ and m^ 2 ) which may have different masses and thus a non- vanishing mixing 
angle occurs. The super partners of the n q light quarks are assumed to have degenerate masses 
(m,q) and vanishing mixing angle. A generalization to a non-degenerate spectrum is possible in 
a straightforward way from the formalism for the top squark sector which is discussed in detail 
in the following. The gluino mass is denoted by rrig. 

Most of the formulae which we are going to present in the following can already be found in 
Ref. [5]. For completeness we repeat the most important ones here and extend them to three 
loops. Unless stated otherwise all parameters in the following derivation are DR quantities which 
depend on the renormalization scale /i. For the sake of compactness the latter is omitted. Bare 
quantities are marked by a superscript "(0)". 

It is common to denote the left- and right-handed components of the top squark by ti and t^, 
respectively. The corresponding mass matrix is given by 



with X t = A t — /^susy cot (3. A t is the soft SUSY breaking tri-linear coupling, and Mjy and Mq 
are the soft SUSY breaking masses. With the help of the unitary transformation 







(2) 



it is possible to diagonalize .M- 




(3) 



where the eigenvalues are the masses of the eigenstates t\ and t2- They read 




(4) 
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The unitary transformation can be parametrized by the mixing angle 



cos,, -an* (5) 
y sm&t cosVt J 

with 

sin (2^) = i g 2 • ( 6 ) 

m~ — m~ 

The renormalization constants connected to the top squark are extracted from the top squark 
propagator. At tree-level it is a diagonal 2x2 matrix which receives non-diagonal entries at 
loop-level. It is convenient to absorb the corresponding counterterms into a renormalization 
constant for the mixing angle which we introduce via 

ef ] = e t + 6o t . (7) 

In order to be able to write down the renormalized top squark propagator we define the renor- 
malization constants as follows: The wave function renormalization constant defined through 

(?H /2 (t) 

~ 1/2 

can be parametrized by a universal factor Z 2 and the renormalization constant for the mixing 
angle 

r i/2 _ ^1/2 ( cos 59 t sin 59 t 
% cos 59 t 



* z > —smd9 f cos 69t 1 



This equation follows from Eq. (J2]) and (t^\i^) T = Z^ 2 (tL,tji) T . Furthermore, the renormal- 
ized mass matrix can be parametrized as follows 



(™f) 2 \^(m\ 1 Z mil m\ 2 Z mi2 



{mf^f ) V m 2l Z m 2 i m22 Z m 22 



= M , (10) 



where we require that the off-diagonal elements in the renormalized mass matrix vanish. As a 
consequence, the counterterm 59t takes care of the divergences in the self-energy contribution 
where a t\ transforms into a t<i or vice versa. This can be seen in the explicit formulae given 
below. The diagonal elements of Eq. (|10p can be identified with the renormalization of the 
masses 

O t ~ 0) ) 2 = rn&Zw = m\Z m . (11) 

In order to formulate the renormalization conditions it is convenient to consider the renormalized 
inverse top squark propagator given by 

iS"V) =P 2 (4 /2 )V 2 - W - nP 2 )] Zf (12) 



4 



where 



S n (p 2 ) S 12 (p 2 ) 
S 2 i(p 2 ) S 22 (p 2 ) 



(13) 



stands for the matrix of the squark self energy. In the DR scheme the renormalization conditions 
read 



'j 



0, 



(14) 



pp 



where "pp" stands for the "pole part". 



In order to obtain explicit formulae for the evaluation of the renormalization constants it is 
convenient to define perturbative expansions of the quantities entering Eq. (fill) . Up to three- 
loop order we have 



* = i + (£K> + (£)'«f + (£)V> + o ( 



a , 



\4irJ jj 



AttJ v 



4tt 



,(3) 



+ 0(a 



4 ) 



(15) 



where S {1 5 2} and k € {2,mj ,77i£ }. Inserting these equations into (|12p one can solve 
Eq. (I14p iteratively order- by-order in a s . At one- loop order one gets 



(1) 



(1) 

2 



1,2, 



E£ - <^ (1) ( ml 



'2 



PP 



PP 



(16) 



The terms proportional to p 2 in the first equation of (|16p are used to compute the wave function 
renormalization constant which is independent of all occurring masses. Thus they can be set to 
zero and one obtains 



SZ. 



(i) 



1 

P 2 



pp 



p 



(17) 



pp 



Once 5Z^ is known Eq. f)16|) is used to obtain SZ^~ keeping the mass dependence in 
(see below for more details). The second equation of (|16p is used to obtain the renormalization 
constant of the mixing angle via 



86 



(i) 



J 12 



2 2 

— m~ 



(18) 



pp 



Proceeding to two loops we obtain the equations 



5 



2^ + 59^ (ml -ml 



0, 



1,2. 



pp 



(19) 



59\ ( m i ~ m 



59?5Z^ 



m a _ 2 \ _ rf^a) m ? + 

*1 *2 / 1 m ti *1 1 m t 2 



I'll 



0, 



i']> 



~(2) (2) (2) 

which are solved for Z\ , oZ^. and 59; using the same strategy as at one-loop level. 
Similarly, at three-loop order we have 

2 



) + m ? _ J^W m | - e 



(1) , 



+ 



5 (1) 



}( 2 ), 



5^(m|-m|)-2^S^-2S 12 



(i) v (i) _ 2X3g) 
( 2 ) v (!) 



+ 5Z^ ( E « - 54?) m|) - *Z?^ZW m| + *Z®E 



2^Eg 



0, i = 1,2, 



pp 



(20) 



(21) 



59\ 



(i) 



m t! ""tl 



- <5zf(m| - m|) - tfZ^m? + 4 



3 ( 2 ) 



JzJ 1 ) (m| - ' 



(i) v (i) 



? 1 - SZW m 2 ~ + tfzM ?77 - + S (1) - E (1) 



2 2 
777, r — m_ 



777 r — 771. 

V tl *2 



2 ) 



21591' 



(i)\ v(i) 



J 12 



0. 



pp 



(22) 



Sample diagrams contributing to En up to three loops can be found in Fig. [TJ the contributions 
to E12 and E22 look very similar. Once the quantities En, E12 and E22 are known to three-loop 
order it is possible to extract the renormalization constants for the squark wave function and 
mass and the mixing angle from Eqs. f)21 1) and ()22|) . 

As compared to the corresponding self-energy contributions for fermions or gauge bosons, which 
after proper projection only lead to logarithmically divergent integrals, the quantities in the 
above equations have mass dimension two. As a consequence the renormalization constants of 
the squark masses and the mixing angles depend on the occurring masses, even in a minimal 
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Figure 1: Sample diagrams contributing to En at one, two and three loops. The symbols t, ti, 
g, g and e denote top quarks, top squarks, gluons, gluinos, and e scalars, respectively. 

subtraction scheme like DR. At three-loop order an exact evaluation of the corresponding 
integrals is not possible. It is nevertheless possible to reconstruct the complete dependence 
on the occurring masses using repeated asymptotic expansions and in addition some knowledge 
about the structure of the final result. The latter can be induced from the known results at one- 
and two-loop order. Besides the polynomial dependence inverse powers of first (second) order 
in ~ m i 2 occur in the two-loop contributions to Z m . (59t). Thus we expect that in 5Zm]_ at 

most — m~ ) 2 and in 59^ at most l/(m~ — m~ ) 3 terms appear. Asymptotic expansion 

leads to results where these denominators are expanded in a geometric series. If sufficient terms 
are evaluated it is straightforward to properly reconstruct the inverse mass differences. 

Using asymptotic expansion for several different hierarchies it is possible to check that the 
final result is independent of the actual choice. In our calculation we have chosen the external 
momentum as the largest scale in order to avoid infrared divergences^ and the e-scalar mass as 
the smallest. As far as the squark masses, the gluino and the top quark mass is concerned any 

1 Note that there are still massless gluons and light quarks in the theory. 
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hierarchy can be chosen. We decided to consider the three choices 

q 2 » m| 2 » mj » m~ » to 2 , » m 2 » m 2 , 
<7 2 > mj 3> m~ » m| 3> m~ » m 2 » ?n 2 , 

<? 2 > m| > m~ » m| > m\ > m| > m\ . (23) 

We have checked that in all cases we obtain the same results for Z m . and 50t. Note that in the 
last hierarchy the top quark mass is even larger than the corresponding squark mass which is 
allowed since the mass dependence in the DR counterterms has no physical meaning. 

In each hierarchy of Eq. ([23]) six mass ratios appear. Some of the expansions are simple and can 
be truncated after a few terms. E.g., all terms with inverse contributions in q 2 can immediately 
be set to zero. Similarly, all mass ratios where one has a top squark mass in the denominator and 
mt, rrig or niq in the numerator only low-order expansion terms appear in the final result. This 
has been checked by increasing the expansion depth and verifying that the higher order terms are 
zero. Due to the occurrence of l/(m~ — m~ ) terms in the exact result several terms in m f - /m.- 
have to be kept in the expressions for the self energies in order to be able to reconstruct the 
geometric series. In practice we compute terms up to (mt 1 / m t 2 ) 8 an< ^ check that after including 
two more powers in the top squark mass ratio the final result does not change. 

At this point some comments on the treatment of the e scalar mass, m e , are in order. In 
practice there are two renormalization schemes for m t which are frequently used, the DR and 
on-shell scheme. In the latter one requires that the renormalized mass vanishes in each order in 
perturbation theory whereas in the DR prescription only the pole parts are subtracted by the 
renormalization constant. We will present our results in a first step for DR e scalar masses and 
afterwards discuss the difference to the on-shell scheme. 

In the DR scheme it is important to keep m t different from zero since the renormalization 
group equations for the squark masses and m e are coupled. A non-vanishing e-scalar mass in 
intermediate steps is also required for the computation of the anomalous dimensions in the DR 
scheme [32] (see below) which was constructed in order to disentangle the running of m e from 
the one of the squark parameters. 

After the calculation of the bare self energies we renormalize all occurring parameters in the DR 
scheme. For our three-loop calculation we need the counterterms for a s , mt, m g , ., 6t and m e 
to two-loop order and the one for rriq to one-loop approximation. Furthermore, also the QCD 
gauge parameter has to be renormalized to two loops since it appears in the results for the wave 
function anomalous dimensions. All relevant counterterms can be found in the Mathematica file 
provided together with Ref. [1] and in Ref. [5]. The two- loop corrections for the e-scalar mass 
renormalization is provided in Ref. [33] . 

For the calculation of the three-loop integrals we make use of several computer programs which 
work hand-in-hand in order to reduce the manual interaction to a minimum. All Feynman 
diagrams are generated with the program qgraf [34j . The generated files are manipulated by a 
perl program pQ, which implements the prescriptions of Ref. |35| . in order to obtain the correct 
prefactors due to the Majorana character of the gluino. Afterwards the output is transformed to 
FORM [36] notation with the help of q2e and exp [37 1 138 ] . exp furthermore applies the asymptotic 
expansion (see, e.g., Ref. [39]) in the hierarchies specified in Eq. (|23p . As a result only one-scale 
integrals up to three loops appear which can be evaluated with the packages MINCER [40] and 
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MATAD [41] . Let us mention that we implemented the DRED Feynman rules for SUSY QCD as 
given in Ref . [MJ [Q . 

Once the renormalization constants are available we compute the corresponding anomalous 
dimension with the help of 

dX 



1X X ' 

where the quantity X is either a mass parameter or the mixing angle 



(24) 



X G <{ m~ , m~ , to? m- g , m t , m 2 , 9 t 



(25) 



In practice the derivation in Eq. (|24p is taken after exploiting the relation between the bare 
and the renormalized quantity. Since bare parameters do not depend on [i the derivative acts 
only on the renormalization constant. In the case of the top quark and the gluino the latter are 
mass independent and thus the derivative w.r.t. /i can be rewritten into a derivative w.r.t. a s . 
For the other parameters, however, one has to take into account the mass dependence of the 
Z factors. Let us as an example consider the anomalous dimension of m^. which leads to the 
following chain of equation^! 



7n 



/'" 



Z n 



/'■" 



d 
di? 

dZ 



Z 



da 



. da s x dZ mi . dX 
d^2 + V dX l^fl 



x 



d 
dX 



(26) 



where @{pt s ) is the anomalous dimension of the strong coupling and X runs over the parameters 
listed in Eq. ([25]). 

In the next Section we provide results for various anomalous dimensions. For this purpose it is 
convenient to introduce the following expansion 



lx 



a s \ ^ f ot s 

7T ' V 7T 

n>0 



Tx 



(27) 



3 Results 

In a first step we have computed the three-loop corrections to the squark wave function renor- 
malization constant Z2 (which is mass independent). In the following we present results for the 
anomalous dimensions 7^- , 7 mj , j m . and r ye t up to three-loop order which all have a non-trivial 
mass dependence. The corresponding results for the renormalization constants can be found in 
Mathematica format on the internet page |33j . 

2 In the subscript for the anomalous dimensions we write rtif 1 instead of m| , etc.. 



9 



At one-loop order we obtain the following results 



=Cf 



„(°) 



1 



m~- 



+ - (l - c 4 t) ( m| - m| j + - m t s 2 t 



S2i 



m~ — m\ 



^ + 4 



(28) 
(29) 



where the abbreviations c n t = cos(n6t) and s n t = sin(n^j) have been introduced and Cp = 
(Ng — l)/(2Nc) is the Casimir operator of the fundamental representation of SU(Nc)- In 
Eq. (|28|) we have given the result for 7 mf . The one for y m ~ is obtained by interchanging 
and mi 2 and replacing 9t by —Qf 

The two-loop coefficients read 



■ m„ m t S2t 



1 



+ C f\ -~rnj--(l-c u ) (ml 



Cf Tf in q 
1 



1 



2 2 \ x 2 i 

m r |--fflj + m§ m t s 2 t 



1 2 , " 

-mi 1 



+ - m g + m ? - + — (1 - c 4i ) ( m fi - \ + - m t -m g m t s 2t 



1 



1 , o 1 



^ + ^ mj + ^ (9 - c 4t ) m| + ^ (7 + c u ) m| - ^ - m 5 m t s 2 t 



(30) 



m t 1 

9 9" C2f — - C2t S2t 

m~ — mf 8 

tl 12 



+ n t 



+ CaC F {— C 2t S 2t - - 2 ^ 2 C 2t } + C\ 

16 2 mf — 

^ 11*2^ 



m- g m t 

9 9- C2i 

mf — ra- 
ti t2 



C2t S2t 



1 



TO- — TO, 
tl t2 



2~ °2t — o C 2* s 2t 



(31) 



where Cyi is the Casimir operator of the adjoint representation of SU(Nq) and Tp = 1/2 the 
trace normalization, nt counts the top squark flavours and n q counts the mass-degenerate squark 
flavours and at the same time the massless quarks. In practice we have nt = 1 and n q = 5, 
however, it is nevertheless convenient to keep the labels arbitrary. 



Let us now come to the three-loop results. The anomalous dimensions for the top squark masses 
are given by 



2 (2) ^3 n 2,1 2 ^ 1 . \ ( 2 

m tClrn\ x = C F <3m~g + - m t - -m g m t S2t + — (1 - c 4 t) [m fi 



m 



t-2 



45 



15 



9 



+ C a Cf < —m + — m 5 + -to, - - m~ q m t s 2 t + 777 (1 - c 4t ) m 



32 



+ C F Ca 



16 



4 9 ' 8 4 8 -s— ■ 64 



21 2 3 2 9 3 

— m- to. -| — m?, m+ s 2 t 

8 9 8 * 8 9 64 



ti 



m 



t> 



(1 - c it ) (m\ - m|) 
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C F Tjlnl 
3 



lm 2 e -^m 2 g + jm 2 h -m 2 t +^ m g m t s 2t - ^ {13 - c 4t ) (m? - m| 



3 
4 



1 

32 



+ < 



+ n 9 n t 



2 2 2 

■ m t m= H — 

« 2 9 4 9 4 



- m? + - m§ m< s 2t - — (1 - c 4t ) (m? - m|) 



- m £ - 3m § + - m ? + - m fi - - m t + - m g m t s 2t 



\i7-ou) (ml -ml) 
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16 

- -m g m t s 2t C3 + (m^ - m 



3 

| "»e - — m| + ^ m| - - m z + 3raj m t s 2 i + 9m| C3 + ^ m 2 t ( 3 



1 



13 3 

»J V8 C4 '"2 + 16 C3 " 16 C4iCs 



+ n„ 



m * ~ t m ^ + 1 m ^ ~ m * + 3 m ^ mt S2i + 9 m ^ 3 + \ m * ^ ~ \ m§ mt S2t ^ 



m| - m| 



1 1 3 3 

C4t - - + — C3 - 777 C4t C3 



8 16 



16 



g m t ~ 2 m e - y m l - g m 9 m < s 2t - 9 m? C3 - g m ? C3 



H — rrifl m+ S9f C3 + fwif — m~ , 

2 9 V *i V V 64 64 



11 3 3 

C4t - — C3 + rrz c 4t C3 



16 



16 



+ n t 



2 3 2 15 2 3 

— m c rrir 

2 e 8 *i 



3 

2 m z - ^ m z - m\ - ^ m g m t s 2t - 9 m~ C3 - ^ C3 + ^ rag "it s 2 t C3 



2 9 

+ I mr — mr , 



61 1 3 3 

C4t - — C3 + rrr C4t C3 
16 16 



(32) 



where £3 is Riemann's zeta function with the value £3 = 1.2020569 .... The three-loop expression 
for jg t reads 



(2) 



C F Tf{n 



m g mt 3 

— 2 T 7 c 2* 

mf — m~ 4 



1 

16 



C2i S 2 t 



+ n 2 



m g mt 3 

— 2 r 7 C2 * 

mf — m~ 4 

*1 *2 



1 

16 



C2i S2t 



+ n q n t 



+ Cl 



m„ mt 3 



— 2 T n C2t ~o C 2* S 2t 

mi — mi 2 8 



C2t «2i 



t 2 

m g mt 3 
m~ — m~ 2 

El *2 



™ C 2t s 2i " 2 

62 mi — m 



m g mt 9 



ta 



2 I 9 3 

+ C A Cp \ _ 2 _2 C 2* — 00 C 2* s 2t 



— m~ 

tl *2 



32 



11 



+ C 2 F T f \n q 



m 



m-gmt I 9 \ / 3 1 

2" c 2t ( 3 - - C3 ) + c 2 t s 2t ( - C3 - 7 



m 



rrig m t 




m~ — m, 

tl 12 



+ nt 



mg m t 



2 2 i 

mi — mi \ I 



ti 



t 2 



a 3 ^ ( i 3 ^ 

C3-g l+C 2tS2t (^-gC3 



2 ■ C 2 i 3 - - (3 + C 2t S 2 t o C3 - T 

mi — mi \ 2 \8 4 



mjmi / 9 3\ / 1 3 

2" C 2t ( - C3 - 7 ) + C2t S2t [ — ~ ~ C3 



32 



(33) 



At that point a brief comment on degenerate top squarks is in order. In the expressions for j n 
the limit m 



can be taken naively. Furthermore one has to nullify the mixing angle. 
The quantity jg t is not defined in the mass-degenerate case which is reflected by the fact that 
the limit m ( - 2 — > mr does not exist in Eqs. ()29|) . (|31|) and ()33p . 

In order to compare with the results in the literature we have to transform our results to the 
anomalous dimensions for the quantities Mq , and A t as given in Eq. (jTJ) . This is conveniently 
achieved with the help of Eq. ([3]) which is differentiated w.r.t. fi 2 . The resulting equations are 
then solved for the 7m- , IMq an d JA t - We have compared the resulting one-, two- and three- 
loop expressions with the results in the literature [241 [26l [22] and found complete agreement. 
Note that the method used in Ref. [22] is based on a relation of the anomalous dimensions to 
an all-order expression in the so-called NSVZ scheme [15] whereas in this work a diagrammatic 
approach has been used to evaluate the three-loop corrections. We refrain from providing explicit 
results for JMq an d 7M^ which, however, can be found in the Mathematica file [33]. Note that 
we have JMq — 7a% which is expected since electroweak effects are neglected [29]. This serves 
as a welcome check for our calculation. The result for jA t is proportional to the gluino mass 
and is thus quite compact. Up to three- loop order it is given by 




+ 



a s \ 3 



7T 



-C A C F -C 2 F -C F (n q + n t ) T f 



lc 2 A C F -\c A C F + *-C F --C F (n q + n t ) 2 T 2 



+ C A C F (n q + n t ) T f Q - ^ Cs) + Cf (n q + n t ) Tj ("3 + ^ Cs) j 



(34) 



For completeness let us also provide the result for mass-degenerate squarks which is given by 
,2_(o) =n „ m 2 (35) 



m-qT^ = C F m g , 
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r 2 -m 2 
°F 7 m 5 



C F T f n, 



1 o ^ 9 2 

- m e + - m g + m q 



+ n t 



3 

2^9 



mi H — m\ H — m~ H — ?nj — m 



, (36) 



12 



21 



45 



m 



hml = ci*m\ -c A cii^ m i + ^ m n + c\c F r^ m i + ^ 



16 



32 



15 2 

— fn~ 
4 9 



+ C F Tf{nl 



3 o 3 o ^2 

— mi m~ H — m- 

8 6 2 9 4 9 



3 

- — 3 m~ + - m~ H — m| H — m~ — -ra, 
4 e 9 '4 9 'g ti 'g t 2 4 { 




3 9 3 9 3 9 *^2 2 

— m, m h -\ — m r H — mr — -m f 

e 2 fl 8 *! 8 * 2 4* 



2 2 2a 

- m e + — m~ q + 9 C3 



3 o 1^ 2 2 



1^ 2 2 > 

— m t + 9mgC3 
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(37) 



One observes that all terms which do not involve nt can be obtained from j m . by setting 



■m 



m 



ti ' 



= and 6 t = 0. 



When applying the anomalous dimensions derived in this paper one has to consider the combined 
set of differential equations of all DR parameters appearing on the r.h.s. of the above results. 
This concerns in particular the unphysical e-scalar mass which means that although m e is set 
to zero at one scale it is different from zero once this scale is changed. A way out from this 
situation is to renormalize the e scalar mass on-shell. We have computed the resulting anomalous 
dimensions and provide the results in Ref. [33]. Alternatively one could shift the squark masses 
by a finite term which is chosen such that the e scalar decouples from the system of differential 
equations. The resulting renormalization scheme is called DR' scheme and has been suggested 
in Ref. [32]. In our approximation the finite shift is needed up to two loops which is given 
by [321 M 



2 9 

mj mj 



1 



7T 2 



C F mt + 



C F m e 



1 13 

- T f (n q + n t ) + -C F - -C A 



(38) 



where f = t or f = We have checked that after inserting this shift in j m . , -y m - and 7 m _ the 
parameter m e drops out from the resulting anomalous dimension. Again we refrain from listing 
explicit results, however, provide the analytic expressions in [33]. 



All results presented above can be found in Mathematica format on the webpage [33]. In 
addition we provide the results for the anomalous dimensions 7m- j and 7^ t and the renor- 
malization constants for the squark masses and the mixing angle in the top squark system. The 



Of course, Tf is not altered. 
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1 loop 


2 loops 


3 loops 


m fl (GeV) 


1425 


1416 


1378 


m h (GeV) 


1677 


1670 


1632 




0.658 


0.659 


0.656 


mg (GeV) 


1580 


1573 


1535 



Table 1: Numerical values for the DR.' parameters for \i = Mz using the numbers in Eq. (f39j) as 
input and solving the system of differential equations with one-, two- or three-loop anomalous 
dimensions in the squark sector. 

Mathematica file contains furthermore the result for 7 m , , Jm i2 and y m ~ for on-shell e scalar 
masses and in the DR scheme. 

Let us finally perform a simplified analysis in order to exemplify the numerical impact of the 
three-loop corrections. In our example we fix the following values of the DR parameters at the 
scale 11 = fi G = 10 16 GeV 

m h = 400 GeV , m t = 67 GeV , 9 t = 0.1 , a s = 0.0425 , 

TO£ 2 = nig = rrig = 600 GeV , (39) 

and use the anomalous dimensions obtained in this paper and in Ref. [1] to compute the cor- 
responding values for fj, = Mz- Since our aim is to study the numerical importance of the 
three-loop anomalous dimensions in the squark sector we neglect all threshold effects. Further- 
more, we use for the running of a s , m g and mt always the three- loop approximation whereas in 
the case of the squark masses and Q% the loop-order is varied from one to three. 

The values of mt = mt(fiG) an d ot s = a s (/ic?) in Eq. (|39|) are chosen such that three-loop running 
leads to mt(Mz) = 170 GeV and a s (Mz) = 0.118. The results for mi , rnx, 9t and m^ at the 
scale (j, = Mz can be found in Tab. [TJ 

We observe a small change in the mixing angle by about 0.4%. As far as the squark masses are 
concerned one observes a moderate shift of a few GeV when going from one to two loops. After 
switching on the three-loop terms, however, the squark masses are decreased by about 40 GeV 
which is approximately an order of magnitude larger than the two-loop corrections. Nevertheless 
it corresponds to a shift in the masses of about 3% which is a reasonable amount for a three-loop 
SUSY QCD term. Our observation coincides with the findings of Ref. [22] where also relatively 
large three-loop corrections for the squarks have been identified. 

4 Conclusions 

In this paper the renormalization constants for the squarks and the corresponding mixing angle 
have been computed to three-loop order within super symmetric QCD. Thus, all anomalous 
dimensions of the physical parameters are now available to order a$ and can thus be used to 
relate their mass values at the GUT and electroweak or TeV scale. 

Our calculation has been performed using dimensional reduction for the regularization of the 
divergent loop integrals which is realized with the help of massive e scalars. As far as the 
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renormalization of the e scalar mass is concerned we have evaluated our results for three different 
schemes: DR, OR' and on-shell. Our results agree with Ref. [22] which supports the consistency of 
DRED with SUSY QCD since in Ref. [22 j the results have been obtained without a diagrammatic 
calculation. 

A simplified numerical analysis shows that the three-loop corrections to the squark masses are 
numerically important (see also [22]) and thus should be included in the spectrum generators 
which incorporate the running from the GUT to the electroweak scale. 

All renormalization constants and anomalous dimensions computed in this paper can be down- 
loaded from the URL [33] in Mathematica format. 
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